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摘 要 : 针对 当前 大 部 分 光伏 (photovoltaic，PV) 模 型 参数 辨识 算法 均 存 在 准确 性 低 和 可 靠 性 差 等 问题 ， 提 出 了 一 种 
采用 改进 型 共生 生物 搜索 算法 (symbiotic organisms search，SOS) 的 光伏 组 件 模 型 参数 辨识 方法 。 首 先 ， 为 提高 标准 
SOS 算法 的 寻 优 性 能 ， 提 出 了 新 的 改进 型 SOS 算法 ， 记 作 ImSOS 算法 。 该 算法 在 标准 SOS 算法 的 生物 种 群 初 始 化 
阶段 采用 了 准 反射 学 习 机 制 ; 在 互利 共生 搜索 阶段 采用 了 改进 受益 因子 策略 ; 在 偏 利 共生 搜索 阶段 采用 了 收缩 随机 
数 产 生 因子 区 间 策 略 。 其 次 ， 给 出 了 采用 所 提 ImSOS 算法 求解 基于 实验 测量 电流 一 电压 (LV) 数 据 的 光伏 组 件 模型 
参数 辨识 问题 的 具体 步骤 及 实现 流程 。 最 后 ， 利 用 实际 Sharp ND-R250A5 光伏 组 件 进行 实验 ， 通 过 与 标准 SOS 算 
法 以 及 其 他 七 种 新 疾 智 能 优化 算法 进行 对 比 验证 ,结果 表明 了 所 提 ImSOS 算法 在 光伏 组 件 模 型 参数 辨识 的 有 效 性 和 
优越 性 。 可 见 所 提 ImSOS 算法 为 准确 可 靠 地 办 识 光伏 组 件 模型 参数 提供 了 一 种 新 的 有 效 方法 。 
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Parameter identification of photovoltaic module models using improved 
symbiotic organisms search algorithm 
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Abstract: To solve the disadvantages of the most photovoltaic (PV) models parameter identification algorithms at present, 
which have low accuracy and poor reliability, this paper proposed an improved symbiotic organisms search (SOS) algorithm 
for parameter identification of PV module models. First, to enhance the performance of original SOS, a novel improved 
SOS algorithm, named as ImSOS, was proposed. In ImSOS, a quasi-reflection-based learning (ORBL) scheme was 
employed in the population initialization step of original SOS. Moreover, the strategy of the modifications of benefit factors 
was used in the mutualism phase of SOS. A strategy of narrowing the search range of randomly generated coefficients was 
adopted in the commensalism phase of SOS. And then, the procedures and flowchart of employing the proposed ImSOS for 
solving the PV module models parameter identification problem based on experimental current versus voltage (I-V) data of 
a real PV module was detailed. Finally, the proposed ImSOS was demonstrated on the parameter identification of different 
PV module models of the Sharp ND-R250A5 PV module. Experimental results and comparisons with original SOS and the 
other seven novel intelligent optimization algorithms implied the effectiveness and superiority of the proposed ImSOS. 
Therefore, the proposed ImSOS becomes a new effective method to accurately and reliably identify PV module models 
parameters. 
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0 ”引言 大 发 电 系统 的 核心 部 分 ， 同 时 也 是 最 重要 的 部 分 ， 光 伏 组 件 
B 可 以 实现 将 太阳 能 直接 转换 为 电能 Bl]。 因 此 ， 对 光伏 组 件 进 
能 源 一 直 是 人 类 经 济 与 社会 发 展 的 基础 , 特别 是 进入 21 行 准确 地 建 模 与 仿真 对 光伏 发 电 系统 的 设计 计算 以 及 性 能 优 
世纪 以 来 ， 我 国 经 济 水 平 的 迅速 发 展 使 得 对 电力 能 源 的 需求 。 ”化 显得 尤为 重要 。 
大 幅 增 长 。 因 此 ， 持 续 快 速 增长 的 电力 需求 与 传统 化 石 能 源 光伏 组 件 的 建 模 有 三 个 主要 过 程 。 首 先 选 择 合 适 的 等 充 
的 枯竭 之 间 的 矛盾 日 益 凸 显 ， 而 且 出 于 环境 保护 和 可 持续 发 。 电路 模型 ， 其 次 是 建立 其 数学 模型 表达 式 ， 最 后 确定 模型 中 
展 的 要 求 , 这 极 大 地 推动 了 可 再 生 能 源 资 源 的 开发 与 利用 口 ， 的 未 知 参数 。 由 于 光伏 组 件 一 般 是 由 多 个 光伏 电池 串 并 联 构 
特别 是 其 中 的 太阳 能 资源 。 鉴 于 其 绿色 环保 、 可 再 生 及 取 2 成 的 ， 而 实际 中 单 二 极 管 模型 和 双 二 极 管 模型 是 光伏 电池 建 
不 尽 等 诸多 大 势 ,近年 来 通过 光 大 发 电 系 统 来 利用 太阳 能 资 。 模 中 常 采 用 的 两 种 等 效 电路 模型 息 ， 所 以 这 两 种 光伏 电池 模 
源 的 方式 已 在 全 球 范围 内 引起 了 极 大 关注 外 。 然 而 ， 作 为 光 ”型 也 是 建立 光伏 组 件 模型 的 基础 ， 而 且 描 述 这 两 种 模型 的 电 
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流 一 电压 (LV) 特 性 方程 是 复杂 的 隐 式 非 线 性 方程 ， 参 数 多 且 受到 以 上 启发 ， 针 对 已 有 光伏 模型 参数 辨识 研究 中 所 存 
求解 难度 很 大 喇 。 因 此 ， 考 虑 到 上 述 情形 ， 精 确 地 进行 光伏 ”在 的 问题 , 本 文 提 出 了 一 种 采用 改进 型 SOS 算法 的 光伏 组 件 
组 件 模型 中 未 知 参数 的 辨识 是 进行 太阳 能 光伏 系统 设计 与 仿 ” 模型 参数 辨识 方法 ,首先 , 为 提高 标准 SOS 算法 的 寻 优 性 能 ， 
真 计 算 、 性 能 优化 与 评估 、 故 障 检测 与 诊断 、 最 大 功率 点 跟 ”提出 了 新 的 改进 型 SOS 算法 ， 记 做 ImSOS 算法 。 该 算法 在 
踪 控 制 以 及 发 电 输 出 功率 预测 的 重要 基础 和 关键 所 在 ， 这 是 。 标准 SOS 算法 的 生物 种 群 初始 化 阶段 采用 准 反射 学 习 机 制 
一 项 至 关 重 要 而 又 非常 具有 挑战 的 工作 ， 同 时 也 是 当前 的 热 以 提高 算法 的 求解 精度 和 加 快 算法 的 收敛 速度 ， 在 互利 共生 
点 研究 方向 之 一 [6 71。 搜索 阶段 采用 改进 受益 因子 策略 以 使 每 个 生物 在 它们 相互 作 
近年 来 ， 很 多 学 者 提出 各 种 方法 来 解决 光伏 模型 参数 辩 ”用 中 获得 相等 的 部 分 益处 来 提高 算法 的 稳定 性 ;在 偏 利 共生 
识 问题 ， 主 要 分 为 解析 法 、 数 值 计算 方法 和 元 启发 式 优化 方 搜索 阶段 采用 收缩 随机 数 产 生 因 子 区 间 策 略 ， 把 搜索 范围 控 
法 三 大 类 。 文 献 [4] 提 出 了 采用 解析 法 来 提取 基于 实测 LV 数 制 在 合理 区 间 来 加 快 算法 的 收敛 速度 的 同时 也 保证 其 求解 精 
据 的 太阳 能 电池 的 单 二 极 管 模型 和 双 二 极 管 模型 参数 ， 文 献 ” 度 。 其 次 ， 给 出 了 采用 所 提 ImSOS 算法 求解 基于 实验 测量 
[8] 提 出 了 先 利用 泰勒 级 数 展开 式 将 隐 式 单 二 极 管 模型 建 模 。 LV 数据 的 光伏 组 件 模型 参数 辨识 问题 的 具体 步骤 及 实现 流 
成 显 式 解析 模型 ， 进 而 对 太阳 能 电池 的 两 种 改进 五 参数 模型 — 程 。 最 后 ， 利 用 实际 Sharp ND-R250A5 光伏 组 件 进行 实验 ， 
进行 参数 提取 的 方法 。 尽 管 解析 法 能 快速 求解 ， 但 其 求解 适 通过 与 标准 SOS 算法 以 及 其 他 7 种 新 颖 智能 优化 算法 进行 对 
应 性 差 ， 特 别 是 该 方法 中 的 近似 处 理会 大 大 降低 求解 的 准确 比 验证 , 结果 证 明了 所 提 ImSOS 算法 在 光伏 组 件 模型 参数 辩 
度 。 文 献 [9] 提 出 了 采用 高 斯 - 赛 德 尔 迭 代 法 来 辨识 基于 单 二 识 的 有 效 性 和 优越 性 。 
极 管 建立 的 光伏 组 件 (single-diode PV module, SDPM) 模 型 参 DET. i ei 
数 。 尽 管 数值 计算 方法 能 提供 准确 的 结果 ， 但 结果 的 准确 程 1 “光伏 组 件数 学 模型 
度 依 赖 于 初始 值 的 选择 。 而 且 上 述 两 类 方法 都 对 模型 有 一 定 1.1 光伏 电池 模型 
的 限制 条 件 ( 如 要 求 在 凸 集 上 连续 可 微 等 ), 存在 不 确定 收敛 、 近年 来 ， 不 同 的 等 效 电路 模型 被 相继 提出 来 描述 光伏 电 
计算 量 大 等 缺陷 。 随 着 人 工 智 能 的 莲 勃 发 展 ， 尤 其 是 各 种 新 ” 池 的 电气 特性 ， 其 中 最 实用 而 且 应 用 最 广 的 两 种 模型 分 别 是 
颖 的 元 启发 式 算 法 的 不 断 涌现 ， 其 越 来 越 多 的 应 用 于 光伏 模 。 单 三 极 管 模型 和 双 二 极 管 模型 ， 同 时 这 两 种 光伏 电池 模型 也 
型 的 参数 辨识 问题 中 并 掀起 了 一 股 热潮 ， 如 遗传 算法 (genetic 是 建立 光伏 组 件数 学 模型 的 基础 。 因 此 接 下 来 先 分 别 介绍 描 
algorithms, GA)!” 、 粒 子 群 优化 算法 (particle swarm  ” 述 光伏 电池 特性 的 这 两 种 模型 。 
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(double-diode PV module，DDPM) 模 型 参数 辨识 的 下 
报道 。 因 此, 探索 为 光伏 组 件 模型 参数 辨识 问题 


作 主 要 针对 光伏 电池 模型 参数 辨识 开 


在 某 个 算法 能 
参数 辨识 问题 
展 ， 而 对 光伏 


立 的 光伏 组 件 


高 精度 和 
题 ， 同 时 


可 靠 性 的 元 启发 式 算法 是 一 个 亚 待 解 
也 是 工程 实践 的 热点 。 


由 Cheng fll Prayogo 通过 模拟 自然 界 


究 鲜 见 
寻找 新 颖 的 、 
决 的 关键 性 难 


生生 物 搜索 算法 (symbiotic organisms search,，SOS)[26] 
是 2014 年 


中 不 同 生 物 之 


kt 生 、 偏 利 共 生 和 寄生 三 种 主要 交互 


间 的 互利 
高 自身 适 


应 环境 得 以 生存 的 能 力 而 提出 的 一 种 


式 优 化 算 
群 规模 和 
参 工 作 ， 
[27-29] E ER 


zt 


法 。 该 算法 结构 简洁 ， 操 作 简单 ， 而 
且 优 化 能 力 强 ， 


一 经 提出 便 得 到 了 极 大 的 关注 


作用 方式 以 提 
新 颖 的 元 启发 
且 除 了 生物 种 
进行 繁 元 的 调 


容易 出 现 
多 峰 复 杂 
性 能 。 

对 并 
An AH 
提出 的 一 


优 精度 和 


验证 了 准 
最 优 解 。 


算法 中 来 


习 机 制 (quasi-reflection-based learning, 


早熟 收敛 陷入 局 部 最 优等 问题 ， 特 别 
优化 问题 时 表现 得 更 明显 ， 需 要 进 


学 习 机 制 (opposition-based learning，OBL)B9 最 


铁 卢 大 学 的 Tizhoosh 教授 于 2005 年 
种 新 策略 ， 已 被 广泛 应 用 于 提高 智能 
加 快 其 收敛 速度 。 然 而 最 近 有 学 者 提 


反射 解 相 对 于 对 立 解 而 言 更 有 可 能 靠 


Tj. 与 大 多 数 智能 优化 算法 类 似 , 标准 


QRBL)P 


ESOS 算法 也 
是 在 求解 高 维 
步 改 善 其 优化 


1.1.1 单 二 极 管 模型 


I-V 特性 而 被 广泛 


单 二 极 管 模型 


忆 其 能 简单 和 准确 地 描述 光伏 电池 的 静态 


阻 以 


ifi h 


其 中 
联 电 
Ia 23 


针对 机 器 智能 
优化 算法 的 寻 
学 


出 了 准 反 射 
] 的 概念 而 
近 问 题 的 全 局 
些 智能 优化 


>H 


此 ， 准 反射 学 习 机 制 已 成 功 应 用 于 


提高 它们 的 优化 能 力 2 99, 


图 中 可 以 看 出 ， 该 模型 
电阻 构成 
(Kirchhoff's current law, KCL), 


及 并 联 


为 


Iph 为 光 生 E 


MAJEN. Hesh HENIE h T ERE Rn 


I=L [exp(q(V, +I,R,)/akT)-1] 

其 中 :4 为 二 极 管 反 向 饱和 电流 ;Vt 为 光伏 
Rs 为 串联 
(1.3806503 x 10? J/K) q 为 电子 
T 为 光伏 电池 绝对 温度 。 而 上 且 
(Kirchhoff's voltage law,，KVIL)， 可 知 流 过 并 联 电阻 的 
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极 管 模型 中 输出 
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由 式 (4) 可 知 ， 在 光伏 电池 的 单 二 极 管 模型 中 待 确 
定 的 五 个 未 知 参 数 分 别 为 Jphn、Jsa、a、Rs 和 Rsh。 
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图 1 
Fig. 1 
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1.1.2 双 二 极 管 模型 从 式 (10) 可 以 清楚 地 发 现 ， 当 光伏 组 件 模型 以 基于 光伏 


在 实际 的 光伏 电池 建 模 中 ， 双 二 极 管 模型 是 除了 单 二 极 “ 电池 双 二 极 管 模型 来 建立 时 , 则 DDPM 模型 中 待 确定 的 未 知 
管 模型 之 外 应 用 最 广泛 的 。 尽 管 单 二 极 管 模型 能 对 实际 光伏 ”参数 增加 至 七 个 ， 分 别 为 Dos. la a2. a1. a2. REI Rsh。 
电池 特性 提供 令 人 满意 的 贴近 度 ， 但 双 二 极 管 模型 考虑 了 耗 。 这 些 未 知 的 模型 参数 同样 地 需要 利用 实测 I-V 数据 进行 辨识 ， 
尽 区 中 复合 电流 损失 的 影响 ， 使 模型 更 加 逼真 并 获得 更 高 的 ”而 且 当 模型 参数 增加 时 ， 其 辨识 难度 更 大 。 因 此 ， 准 确 地 辩 
准确 度 。 因 此 ， 双 二 极 管 模型 等 效 电路 如 图 2 所 示 。 由 图 可 识 这 些 未 知 参数 以 确保 实际 光伏 组 件 有 更 好 的 性 能 ， 这 在 太 
知 ， 该 模型 由 光 生 电流 源 、 两 个 并 联 二 极 管 、 串 联 电 阻 以 及 ”” 阳 能 光伏 发 电 系统 的 相关 研究 设计 中 显得 尤为 重要 。 

并 联 电 阻 构成 时。 同样 , 根据 KCL 可 得 出 图 中 光伏 电池 输 2 ”标准 SOS 算法 及 ImSOS 算法 


出 电流 五 为 
和 = 一 人 (5) 2.1 标准 SOS 算法 
其 中 : a 为 流 过 第 一 个 二 极 管 的 电流 ;Za 为 流 过 第 两 个 二 极 SOS 算法 是 2014 年 由 Cheng 等 人 29 提 出 的 一 种 新 颖 的 
EREM. WI HAKI, Ia Dao 分 别 为 元 启发 式 优化 算法 。 该 算法 源 于 对 自然 界 生 态 系统 中 各 生物 
Tula [exp(a(V, * 1,R.)/a&T) -1] (6) ， 体 之 间 通 过 相互 作用 影响 来 共同 生存 繁衍 策略 的 模拟 。 自 然 
In7l,5 [exp(q(V. * LR. )/aXkT) -1] (7) 界 中 各 生物 体 之 间 最 常见 的 共生 关系 分 为 互利 LÆ, dim 利 
其 中 : al 和 a 分 别 为 第 一 个 和 第 二 个 二 极 管 反 向 饱和 电流 : 生 和 寄生 三 种 方式 。 而 SOS 算法 正 是 通过 简单 的 数学 模型 3 


来 
ai 和 a2 分 别 为 第 一 个 和 第 二 个 二 极 管 理想 因子 。 类 似 于 单 二 模拟 上 述 生物 体 之 间 相 互 作 用 的 生存 方式 实现 寻 优 过 程 ， 这 
极 管 模型 ， 光 伏 电 池 的 双 二 极 管 模型 中 输出 电流 和 输出 电压 样 可 以 用 来 指导 求解 优化 问题 。 因 而 标准 SOS 算法 的 四 个 主 
的 关系 可 改写 为 要 步骤 如 下 。 

L-I, — La [exp(q(V, + LR,)/akT)-1 2.1.1 生物 种 群 随机 初始 化 阶段 


La; [exp(a(V, € LR )/akT) - 1]- (8) 首先 在 搜索 空间 范围 内 随机 初始 化 n 个 生物 如 下 : 
(V, * LR,)/R,, 0,=lb+rand (0, 1) * (ub — Ib) (11) 
从 式 (8) 中 可 以 看 出 ,光伏 电池 的 双 二 极 管 模型 中 待 确定 TX OD 生成 个 初始 解 ， 其 中 : 90; 表示 生态 系统 中 
的 未 知 参数 有 七 个 ,分 别 为 Dons Iais Isas ais azs Rs 和 Rs。 第 i(i=1,2,…,) 个 生物 ， 且 每 个 生物 都 是 一 个 d 维 向 量 ; 
—« rand(0,1) 表 示 一 个 向 量 ， 且 向 量 中 的 元 素 都 是 [0，1] 间 的 随 
全 | 本 | 更 | 汪汪 机 数 ;ub 和 Ib 分 别 为 搜索 空间 范围 的 上 限 向 量 和 下 限 向 量 。 
m VY V 28 v 2.1.2 互利 共生 搜索 阶段 
SOS 算法 中 通过 模拟 自然 界 生 态 系 统 中 如 蜜蜂 跟 花灯 之 
- . 间 的 互利 共生 关系 形成 互利 共生 搜索 阶段 。 在 这 个 阶段 ， 假 
图 2 ”光伏 电池 的 双 二 极 管 模型 等 效 电路 图 设 OW 9; 分 别 表示 生态 系统 中 随机 选择 的 第 i 个 和 第 j (0 
Fig.2 Double-diode model of PV cell ži) NÆ, M 08; 和 6) 两 个 生物 通过 互利 共生 搜索 ， 使 得 
1.2 光伏 组 件 模型 各 自 向 最 优 个 体 学 习 ， 分 别 产生 新 解 9inew 和 0jnew 计算 如 
在 实际 应 用 中 ， 单 体 的 太阳 能 电池 不 能 直接 当 作 电源 使 F: 
用 ， 而 是 将 多 个 太阳 能 电池 串联 和 《或 者 ) 并 联 构成 太阳 能 Oisew =0, trand (0,1) * (6, — MV * BE) (12) 
光伏 组 件 以 提高 其 输出 电压 和 (或 者 09 HB. Fe pL 1920), DR] JI, nox 70; rand (0.1) * (6, - MV * BE.) (13) 
光伏 组 件 模型 的 建立 是 以 光伏 电池 模型 为 基础 的 ， 下 面 分 别 其 中 : MV. BF 和 BF 可 分 别 表示 为 
介绍 两 种 光伏 组 件 模型 ， 即 基于 单 二 极 管 建立 的 光伏 组 件 MV =(0:+0;)/2 (14) 
(single-diode PV module, SDPM) 模 型 和 基于 双 二 极 管 建立 的 BF =round (1+rand (0,1)) (15) 
光伏 组 件 (double-diode PV module, DDPM) 模 型 。 BF,=round (1+rand (0,1) (16) 
1.2.1 SDPM 模型 其 中 : njetL2z,niHje£ü Obest 表 示 生 态 系统 中 当前 最 
当 光 伏 电 池 电 气 特性 由 单 三 极 管 模型 来 描述 时 ， 则 根据 ”” 优 个 体 ( 对 应 当前 目标 函数 值 最 小 的 解 向 量 ); MV 为 互利 向 
KCL 可 得 出 SDPM 模型 中 输出 电流 和 输出 电压 的 关系 为 量 ， 表 示 两 个 生物 体 间 的 关系 特征 向 量 ; rand(0,1) 为 [0, 1] 
LIN, — TaN, [exp(a(V, + 7.R, N.N, JaNkT) -1]- 0) 间 的 随机 数 ，BF 和 BF: ESAKA, HT pap 1 
(V € 4R, NJ N,)/ Rs N/N, 或 2， 这 两 个 因子 表示 每 个 生物 从 相互 作用 关系 中 的 受益 程 
其 中 :Ns 和 No 分 别 为 光伏 电池 串联 和 并 联 个 数 。 BE; Oxa - MV * BF) FI (Oxa -MY*BP) 反 映 了 两 个 生物 在 互利 
很 显然 ， 由 于 光伏 组 件 模型 是 基于 光伏 电池 单 二 极 管 模 ”共生 中 各 自 受益 来 提高 它们 的 生存 优势 并 各 自 趋 向 最 优 个 体 


型 而 建立 的 ， 所 以 在 SDPM 模型 中 需要 确定 五 个 未 知 参数 ， ”学 习 以 达 到 适应 环境 的 能 力 。 在 这 个 阶段 ， 6; 和 267 两 个 生 
分 别 为 pis fa. a. Rs Rsn。 而 生产 厂家 的 数据 手册 一 般 不 ” 物 通过 上 述 机 制 产生 的 新 解 能 否 被 接受 取决 于 它们 的 新 解 所 
会 提供 上 述 参 数值 ， 可 利用 实测 LV 数据 对 这 五 个 参数 进行 。 对 应 的 目标 函数 值 是 否 比 相互 作用 前 它们 对 应 的 目标 函数 值 


辨识 。 而 且 精 确 地 提取 这 些 模型 参数 值 能 使 模型 更 好 地 描述 更 好 。 
光伏 组 件 的 实际 电气 特性 。 2.1.3 偏 利 共生 搜索 阶段 
1.2.2 DDPM 模型 SOS 算法 中 通过 模拟 自然 界 生态 系统 中 如 兰花 跟 乔 木 树 
当 光 伏 电 池 电 气 特 性 由 双 二 极 管 模型 来 描述 时 ， 则 根据 之 间 的 偏 利 共生 关系 建立 偏 利 共生 搜索 阶段 。 在 这 个 阶段 ， 
KCL 可 得 出 DDPM 模型 中 输出 电流 和 输出 电压 的 关系 为 生态 系统 中 随机 选择 的 第 7 个 生物 9; 在 与 第 i 个 生物 0. 
Lol N, - LaN, [exp(a(V, + TR N,/N, Ja NT)-1]- 互 作 用 影响 过 程 中 ， 生 物 8; 从 中 受益 ， 而 生物 29) 不 受益 也 
LoN,[exp(a(V. € LR NVN)J/wNkT)-]- — (10) ”不 受害 。 因 而 产生 的 新 解 9inew 计算 如 下 : 


(V+ LR N,/N, )/ Rs, N.N, 6,70, trand (71,1) * (Owes —6;) (17) 
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其 中 : rand(-1,1) 表 示 一 个 向 量 ， 且 向 量 中 的 元 素 都 是 [-1, 1] 般 而 言 ， 基 于 种 群 的 优化 算法 首先 第 一 步 都 需要 一 些 
间 的 随机 数 ， bus 700 表示 两 个 生物 在 偏 利 共 生 中 ， 由 生物 。 初始 解 ， 然 后 在 搜索 优化 过 程 中 不 断 提高 解 的 质量 以 趋向 最 
90; 提供 的 有 益 优势 来 帮助 生物 91 增强 其 生存 优势 以 不 断 趋 。 优 解 。 这 个 搜索 过 程 会 一 直 重 复 ， 直 到 满足 一 些 预 设 的 准则 
向 最 优 个 体 。 同 样 在 这 个 阶段 ， 若 新 解 对 应 的 目标 函数 值 优 ”” 才 会 停止 。 由 于 对 问题 的 解 缺 乏 一 些 先 验 知识 ， 所 以 在 算法 
于 原生 物 个 体 的 目标 函数 值 ， 则 接受 原生 物 个 体 更 新 产生 的 。” ”初始 化 阶段 多 是 通过 随机 产生 种 群 的 初始 解 。 然 而 有 学 者 已 


$ 


新 解 。 经 验证 了 对 立 解 相对 于 随机 产生 的 近似 解 而 言 拥有 更 高 的 概 
2.1.4 寄生 阶段 率 接近 问题 的 全 局 最 优 解 。 最 近 又 有 学 者 提出 了 准 反 射 学 习 
寄生 关系 是 指 两 种 不 同 生物 的 相互 作用 对 一 方 有 益 ， 对 PORS, 而且 验证 了 准 反 射 解 相对 于 对 立 解 而 言 更 有 可 能 


男 一 方 有 害 。 例 如 症 原 虫 利用 它 与 按 蚊 之 间 的 关系 来 在 人 类 ”靠近 问题 的 全 局 最 优 解 。 
宿主 之 间 传 播 ， 当 症 原 虫 在 人 体内 繁衍 繁殖 时 ， 宿 主人 类 就 为 了 更 清楚 地 解释 准 反射 学 习 的 概念 ， 必 须 先 理解 准 反 
会 遭受 症 疾 甚至 最 后 有 可 能 死亡 。 射 学 习 中 用 到 的 准 反射 数 和 准 反射 点 的 概念 。 

SOS 算法 中 通过 模拟 自然 界 上 述 这 种 关系 建立 寄生 阶段 。 定义 1 假设 x 为 实数 且 rE[a, b], W x f oc PC R 


在 这 个 阶段 , 首先 通过 复制 生物 Or 随机 选取 其 中 某 些 维度 ， 示 为 : X" — rand((a - b)/2,x) , 其 中 rand((a +b)/2, x) 表示 在 (a+b)/2 
替换 为 搜索 空间 范围 内 的 随机 取 值 ， 进 而 形成 一 个 新 个 体 ， 与 x 之 间 均 匀 分 布 的 随机 数 。 

叫 寄生 向 量 , 记 做 PV; 然后 把 从 生态 系统 中 随机 选择 的 第 70 定义 2 假设 P=G,%,…, 羽 ) 为 一 个 n 维 向 量 空间 的 点 ， 
汉 i) 个 生物 6; 作为 宿主 ， 计 算 寄 生 向 量 PV 和 宿主 06; 所 对 其 中 xi, x2, Xn 都 属于 实数 xiE[a; b], Vie{l,2,.…,n} 。 则 
应 的 目标 函数 值 并 进行 比较 。 若 寄生 向 量 PV 所 对 应 的 目标 P RERI A Pe-GgExRe.g xD, Hep rrgs rr ARRA 


函数 值 更 优 ， 则 寄生 向 量 PV 杀 死 宿主 占用 其 位 置 进 而 成 为 xi = rand((a, +b)/2,%) (18) 
第 7 个 生物 体 , 否则 宿主 0; 将 对 其 具有 人 免疫 功能 而 存活 被 保 其 中 :vand((a; +b)/2,%) 表 示 在 (ai+bi)/2 与 xi 之 间 均 匀 分 布 的 
留 下 来 .此 外 , 标准 SOS 算法 的 具体 实现 流程 可 参考 文献 [26]， ”随机 点 。 
这 里 不 再 袭 述 。 这 样 了 解 了 准 反射 点 定义 后 ， 那 么 采用 准 反 射 学 习 的 优 
2.2 ImSOS 算法 化 可 定义 如 下 : 

如 何 权衡 算法 的 全 局 探索 能 力 与 局 部 开发 能 力 决 定 了 算 定义 3 ”假设 P=G,%…, 尺 ) 为 一 个 n 维 向 量 空间 的 点 
法 的 寻 优 性 能 。 通过 分 析 标 准 SOS 算法 可 发 现 其 有 三 个 方 画 UBRE P JIRE, SO 为 候选 解 的 目标 函数 值 ， 此 外 根据 


的 不 足 : a) 标 准 SOS 算法 在 生物 种 群 初始 化 阶段 采用 的 是 随 ” 准 反 射 点 的 定义 可 知 已 的 准 反 射 点 为 Pv -Ggxtesx» 。 若 
机 初始 化 种 群 的 方式 ， 该 种 种 群 初始 化 方式 降低 了 算法 的 全 FPES f(P)， 则 表示 Pr E P RA ERHI HERK m 此 时 
局 探索 能 力 ， 最 终 导致 算法 早熟 收敛 而 陷入 局 部 最 优 解 ，b) ”选择 Ped REP: 否则 保持 了 不 变 。 
标准 SOS 算法 在 互利 共生 搜索 阶段 中 受益 因子 的 选取 方式 为 提高 标准 SOS 算法 的 寻 优 性 能 , 提出 将 准 反射 学 习 机 
容易 MEE TE 可 行 区 域 ， 增加 了 最 优生 物体 与 整个 生 ” 制 引 入 到 标准 SOS 算法 的 种 群 初始 化 阶段 ,通过 考虑 初始 解 
态 系统 之 间 的 目标 函数 的 偏差 ， 进 而 降低 了 算法 的 稳定 性 和 的 同时 也 对 其 准 反射 解 进行 评价 ， 择 优 使 用 ， 这 样 采用 准 反 


可 靠 性 ; c) 标 准 SOS 算法 中 偏 利 共生 搜索 阶段 的 随机 数 产 生 ” 射 学 习 机 制 能 够 对 搜索 空间 进行 充分 全 面 的 搜索 ， 使 得 算法 
姑 子 区 间 容 易 扩大 当前 生态 系统 的 搜索 范围 ， 最 终 使 得 算法 ”快速 收敛 的 同时 又 有 更 高 的 概率 找到 接近 全 局 最 优 解 的 候选 
的 收敛 速度 变 慢 。 解 ， 提 高 了 算法 的 求解 精度 。 

为 克服 标准 SOS 算法 的 缺陷 , 提出 了 新 的 改进 型 共生 2.2.2 采用 改进 受益 因子 策略 的 互利 共生 搜索 
物 搜索 算法 , 记 做 ImSOS 算法 。 该 算法 采用 三 个 主要 策略 来 在 互利 共生 关系 里 ， 参 与 相互 作用 的 两 个 不 同 生物 将 会 
改善 标准 SOS 算法 的 优化 性 能 , 很 好 地 平衡 了 算法 的 全 局 探 ” ”从 它们 的 交互 影响 过 程 中 获得 部 分 或 者 全 部 益处 ， 因 此 在 标 
索 与 局 部 开发 能 力 。 首先, 在 标准 SOS 算法 的 生物 种 群 初始 — ESOS 算法 中 通过 受益 因子 随机 选择 1 或 者 2 来 表征 这 个 特 


化 阶段 采 rm mn 性 .显然 , 若 两 个 生物 从 相互 作用 中 各 自 都 获得 了 全 部 益处 ， 

法 的 收敛 速度 然后， 在 互利 共生 搜索 阶段 采用 改进 受益 因 即 两 个 生物 的 受益 因子 BFi-BFo- 2， 则 上 述 式 (12) 和 (13) 方 
子 策略 以 使 每 个 生物 在 它们 相互 作用 中 获得 相等 的 部 分 益处 ， 程 中 的 MV * BF fll MV * BE, 两 部 分 相 比 在 当 两 个 生物 的 受益 
提高 了 算法 稳定 性 ， 最 后 ， 在 偏 利 共 生 搜 索 阶 段 采 用 收缩 随 AF BFi-BF»- 1 时 ,将 会 在 两 个 生物 之 间 的 关系 特征 中 产生 
机 数 产生 因子 区 间 策 略 ， 把 搜索 范围 控制 在 合理 区 间 ， 加 快 ”更 巨大 的 多 样 性 ， 这 样 直接 导致 OMV BR) 和 
了 算法 的 收敛 速度 的 同时 也 保证 了 求解 精度 。 具 体 改进 如 下 。 Oxa - MV * BE) 两 部 分 产生 很 大 的 多 样 性 , 最 终结 果 就 是 产生 
2.2.1 采用 准 反射 学 习 机 制 的 种 群 初始 化 的 新 解 86inew 和 9jnew 有 可 能 都 不 在 搜索 空间 范围 内 。 因 此 ， 

正如 上 述 2.1 节 中 提 到 的 ， 标准 SOS 算法 中 采用 的 是 随 ” 这 会 使 得 基于 标准 差 评价 的 算法 的 稳定 性 大 大 降低 ， 而 且 算 
机 初始 化 生物 种 群 的 方法 ， 然 而 这 种 随机 初始 化 种 群 的 方法 ”法 收敛 到 全 局 最 优 解 的 速度 减 慢 。 对 此 ， 提 出 采用 改进 受益 
降低 了 算法 的 全 局 探索 能 力 , 因而 导致 标准 SOS 算法 的 收敛 。 ”因子 策略 的 互利 共生 搜索 方式 , 因而 改进 受益 因子 BF1 和 BF， 


精度 下 降 最 终 很 容易 陷入 局 部 最 优 。 对 此 ， 提 出 了 采用 准 反 wF: 

射 学 习 机 制 来 增强 标准 SOS. 算法 的 求解 质量 和 收敛 速度 。 BFi-1 (19) 
对 立 学 习 机 制 最 早 由 加 拿 大 滑铁卢 大 学 的 TizhooshPq 于 BF;-1 Q0) 

2005 年 针对 机 器 智能 提出 的 一 种 新 策略 ,已 被 广泛 应 用 于 提 式 (19) 和 (20) 中 两 个 生物 的 受益 因子 BF 和 BF2, 都 固定 

高 智能 优化 算法 的 寻 优 精度 和 加 快 其 收敛 速度 ， 其 基本 思想 取 值 为 1， 表 征 每 个 生物 在 它们 相互 作用 中 获得 相等 的 部 分 

是 一 个 候选 解 的 对 立 解 《或 称 为 反 向 解 ) 比 该 候选 解 本 身 更 益处， 这 样 新 解 都 能 尽 可 能 在 搜索 空间 范围 内 产生 ， 提 高 

可 能 接近 问题 的 最 优 解 。 因 此 ， 在 对 初始 候选 解 进行 评价 ”算法 稳定 性 的 同时 也 加 快 了 算法 收敛 到 全 局 最 优 解 的 速度 。 
的 同时 ， 也 对 其 对 立 解 进行 评价 ， 然 后 择优 使 用 ， 这 样 可 加 2.2.3 采用 收缩 随机 数 产生 因子 区 间 策 略 的 偏 利 共生 搜索 


快 算法 的 收敛 速度 ， 提 高 找到 最 优 解 的 几率 。 在 标准 SOS 算法 的 偏 利 共生 搜索 阶段 , 通过 分 析 新 解 产 
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生 的 计算 式 (17)， 可 以 发 现 rand(-1,1) 这 个 向 量 对 算法 优化 过 


程 的 收敛 速度 有 着 重 
间 都 是 -1~1 间 的 随机 数 ， 这 村 


in 


要 的 控制 作 / 


围 进 


如 下 : 


其 中 : rand(0.5,0.7) 表 示 
0.7] 间 的 随机 数 。 这 样 通过 收缩 随机 数 产 生 因 子 区 间 在 加 快 


习 ， 该 向 量 中 元 素 的 产生 区 
会 导致 当前 生态 系统 的 搜索 范 
步 扩 大 ， 最 终 使 得 标准 SOS 算法 的 收敛 速度 变 慢 。 为 


此 ， 提 出 采用 收缩 随机 数 产 生 因 子 区 间 策 略 的 1 


Oisew =0; trand (0.5,0.7) * (Oen — 0) 


有 利 共生 搜索 


Q1) 


^i, HIR 


中 的 元 素 都 是 [0.5， 


算法 收敛 速度 的 同时 也 保证 了 算法 的 求解 精度 。 
2.2.4 ImSOS 算法 的 具体 实现 


综 上 所 述 , 为 
了 一 种 新 的 基于 


ImSOS 算法 与 标准 SOS 算法 一 样 结构 简单 、 


提高 标准 SOS 算法 的 优化 性 能 ， 本 文 提出 


述 三 种 改善 策略 的 ImSOS 算法 。 所 提 


且 除 了 只 需 设置 常规 控制 参数 如 种 群 规模 和 最 大 迭代 


代数 外 , 无 须 调 整 任 何 算法 内 部 依 | 
究 者 掌握 和 应 用 。 所 提 ImSOS 算法 的 实现 伪 代 码 如 算 


启 采 用 准 反 射 学 习 机 制 的 种 群 初始 化 阶段 */ 


2.fori=s l:n 
3. forj=1:d 


Mi;- (aj + b;)/2 


4 
5. if (Soi; < Mij) 
6. 


算法 1 所 提 ImSOS 算法 的 实现 伪 代 码 
1. 设置 生物 种 群 规模 n、 


WP TS 


易于 理解 和 实 


制 参数 等 诸多 优点 ， 


最 大 连 代 代数 11wax 和 问题 维 数 4d， 在 搜索 范围 [gj, 妃 内 
随机 初始 化 n 个 生物 种 群 ( 也 就 是 初始 解 集 ) ， 记 作 So。 


UMij: 


QRSo0ij = Soij + (Mij — Soi) x rand(0,1) 


/QRSo: 初始 解 集 So 对 应 的 准 反 射 解 集 ; 


7. else 

8. QRSo;; = Mij + (Soi; — Mi)) x rand(0,1) 
9. end if 

10. end for 

11. end for 


搜索 范围 [w, 世 的 中 点 


rand(0,1): 范围 [0,1] 之 间 的 随机 数 


12. 重新 从 合集 {So, QRSo} 中 选择 最 优 的 n 个 解 作为 初始 解 集 So 


13. while (It Its) 


14. fori2l:n 

15. 确定 当前 最 优 目标 函数 值 对 应 的 生物 体 O pest 

启 采 用 改进 受益 因子 策略 的 互利 共生 搜索 阶段 */ 

16. 对 生物 个 体 9,， 从 生态 系统 中 随机 选择 生物 0, (jui) 
17. BF,-1 

18. BF,-1 

19. MV—(8,* 86/2 

20. 6 ,w= 0, *rand(0,1) * (6... — MV * BF) 

21. Oj new = 8; +rand (0,1) * (Oru — MV * BF,) 

22. 将 新 解 9inew 和 91jwew 分 别 与 9; 和 和 86) 进行 解 的 质量 的 


比较 评价 ， 分 别 选 择 较 优 解 保留 


谍 采 用 收缩 随机 数 产 生 因 子 区 间 策 略 的 偏 利 共生 搜索 阶段 */ 

23. 随机 选择 两 个 不 同 生 物 6; 和 9; gD 

24. Onw = O+rang(0.5,0.7) * (Ope —0,) 

25. 将 新 解 2inew 与 8 进行 解 的 质量 的 比较 评价 ， 选 择 较 优 解 保留 

[* FERE S 

26. 复制 生物 9,， 随 机 选取 其 中 某 些 维度 ， 蔡 换 为 搜索 空间 范围 内 
的 随机 取 值 ， 变 异 得 到 寄生 向 量 PV 

27. 生态 系统 中 随机 选择 的 第 10zD 个 生物 9; 作 为 宿主 

28. 对 寄生 向 量 Py 和 宿主 8) 所 对 应 的 目标 函数 值 进 行 比较 ， 
选择 较 优 的 保留 

29. end for 


30. 天 = 天 +1 
31. end while 
32. 输出 结果 并 可 视 化 


3 ”光伏 组 件 模 型 参数 辨识 的 ImSOS 算法 设计 
3.1 参数 辨识 目标 函数 的 构建 


光伏 组 件 模型 参数 辨识 的 


的 就 是 通过 实验 测 的 IV 数 


据 来 为 光伏 组 件 模型 中 待 辨识 的 参数 寻找 最 优 值 ， 从 而 使 得 


模型 能 准确 地 措 


最 优 参 数值 确定 的 模型 来 计 


述 光伏 组 件 的 实际 输 H 


特性 ， 也 就 是 说 根据 
| 算得 到 的 计算 电流 值 与 实验 测量 
电流 值 之 间 的 误差 尽 可 能 小 。 因 此 ， 通 过 重 写 光伏 组 件 的 


SDPM 模型 的 方程 式 (9) 和 DDPM 模型 的 方程 式 (10) 来 定义 误 


差 函 数 分 别 如 下 : 
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eV, 1,0)=1 aN, — I4N, [exp(q(V, + IR, N,/N, )/aN,kT) -1]- 


(V. - LR, N,J N,)/ Rs N,/N, -1 (22) 
0 — [Ig La a, R, Ra] 
e(V, 1,,0)-15,N, — La N, [exp(q(V, - IR, N./N, )/a.N,kT) -1]— 
LaN, [exp(a(V.  1R, N,/N,)/asN,KT) - 1]- Q3) 


(V HLR N JN, )/Rs N.N, - 1 
0 — [Io a, Lat, a, R, Ra] 
t e(Vi, I 96) 为 误差 函数 ， AX 每 组 测量 电压 与 电流 
值 ， 通 过 提取 的 参数 值 确 定 的 模型 来 计算 得 到 的 计算 电 


H. FH 

数据 

流 值 与 实验 测量 电流 值 之 差 ; 

组 件 模型 中 各 待 辩 识 的 未 知 参数 。 而 
学 模型 为 DDPM 模型 时 ， 

此 ， 


6 为 解 向 量 ， 其 中 元 素 为 光伏 
且 当 光伏 组 件 的 数学 模 
SDPM 模型 时 ， 6=[ Iph, Zsa, a, Rs, Ra]: 当 光 伏 组 件 的 数 
0=[ Iph, Tsd1, Tsd?, a1, a2, Rs, Rshļo Øl 
进一步 构建 以 均 方 根 误差 (root mean square error, RMSE) 
伏 组 件 模 型 参数 辨识 问题 的 目标 函数 定义 如 下 : 


1 N 
Owe(0) = [Soeur y 


Q4) 


数 是 


问题 
目标 
3.2 


据 的 
tik 


4 


参数 的 可 行 性 和 
Sharp ND-R250A5 光伏 组 件 的 I-V 数据 分 别 进行 了 该 组 件 的 


基于 


EL EH: 


N 表示 实验 测量 电压 与 电流 数据 对 的 组 数 。 该 目标 函 
一 个 复杂 非 线 性 超越 函数 ， 其 求解 难度 很 大 。 
本 文通 过 把 光伏 组 件 模 型 参数 辨识 问题 转换 为 一 个 优化 
求解 ， 且 在 搜索 范围 内 不 断 调整 模型 参数 解 向 量 9 以 使 
函数 OFrmsel 0) 最 小 化 。 

采用 ImSOS 算法 的 光伏 组 件 模型 参数 辨识 基本 步骤 
本 节 对 采用 所 提 ImSOS 算法 来 辨识 基于 实验 测量 I-V 数 
太阳 能 光伏 组 件 模型 参数 的 实现 过 程 进行 了 详细 介绍 ， 
步骤 以 流程 图 的 形式 给 出 ， 如 图 3 所 示 。 


算 例 分 析 


为 验证 本 文 所 提 ImSOS 算法 辨识 太阳 能 光伏 组 件 模型 
效 性 ， 现 根据 文献 [7] 实 际 测量 得 到 的 


DE 


单 二 极 管 建立 的 模型 和 基于 双 二 极 管 建立 的 模型 参数 状 


Wi 


光伏 


究 ， 即 算 例 1 光伏 组 件 的 SDPM 模型 参数 辨识 和 算 例 2 
组 件 的 DDPM 模型 参数 辨识 研究 ,该 太阳 能 光伏 组 件 由 


60 块 多 晶 硅 太阳 能 电池 (156.5 mm X 156.5 mm) 串 联 构 成 ， 实 


验 中 
的 最 
路 电 
种 当 
光伏 


[26]. 


声 搜索 算法 (improved harmony search, IHS)U691, MFOU?I, 
MVOP3]、 楼 海 鞘 群 算 法 (salp swarm algorithm, SSA)P4!, 


测试 温度 为 59*°C， 太 阳 辐 照度 为 1 040 W/m?。 在 标准 测 
件 下 (太阳 辐 照度 为 1000 W/m?， 温 度 为 25°C)， 该 组 件 
大 功率 点 电流 为 8.10 A, 最 大 功率 点 电压 为 30.90 V, 短 
流 为 8.68 A， 开 路 电压 为 37.60 V。 同 时 ， 本 文采 用 了 八 
前 较 新 颖 的 人 工 智 能 优化 算法 与 所 提 ImSOS 算法 进行 
组 件 模 型 参数 辨识 性 能 对 比 研究 , 分 别 是 标准 SOS 算法 
带 自 适 应 惯性 权 值 粒子 群 优化 算法 (PSO-w)P1I、 改 进 和 


WOAPR1 及 正弦 余弦 算法 (sine cosine algorithm, SCA)U?l, 为 体 
现 对 比 的 公平 性 ， 本 文 各 算 例 中 所 有 参与 比较 的 辨识 算法 均 


在 CPU X Intel® Core™ 15-2415M., EH 2.30 GHz、 内 存 


为 4 
算法 
算 例 
为 2 


GB 的 计算 机 上 采用 MATLAB 编程 实现 , 而 且 上 述 所 有 
在 各 算 例 中 均 具有 相同 的 种 群 规模 和 最 大 迭代 次 数 ， 即 
1 中 上 述 所 有 算法 均 采 用 种 群 规模 为 10, 最 大 迭代 次 数 
000; 算 例 2 中 上 述 所 有 算法 均 采 用 种 群 规模 为 15， 最 


大 和 返 代 次 数 为 3 000。 算 例 中 各 算法 其 他 算法 内 部 依赖 性 控 


制 参 
各 算 


数 设置 见 表 1。 考 虑 到 人 工 智能 优化 算法 的 随机 性 ， 对 
例 中 上 述 算法 均 独 立 运行 30 次 。 
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HEPER n. IC KGRTVIDCIsses REAd: VEIDUNTH SCUST-V 
数据 中 电压 值 VY 和 电流 值 


Y 
考虑 待 辨识 的 参数 变量 情况 下 初始 化 "个 生物 种 群 〔 也 就 是 初 
集 ) ， 记 作 So。 当 光伏 组 件 3 型 时 ， 则 该 模型 待 辩 
参数 解 向 量 6= Ups RR M 模 型 时 ， 则 该 
辨识 的 参数 解 向 量 OmU, JJ. AERE 6 在 搜 
内 随机 产生 。 


采用 准 反射 学 习 机 制 的 种 群 初始 
化 阶段 


按 式 (18) 产 生 初 始 解 集 对 应 的 准 反射 解 集 ， 记 作 QRSo 


y 
根据 式 (24) 的 目标 函数 OFkmse( 8) 来 计算 初始 解 集 So 中 每 个 解 对 应 
的 人 目标 函数 值 。 


重 SÜBERUES UEFA EQ 选择 最 优 的 n 个 解 作 为 初始 
解 集 5o， 即 把 上 步 中 各 目 e EBORE 选择 前 个 目标 
函数 值 对 应 的 解 为 最 优 的 a 个 解 。 


[er 


t+1 


A | 采用 收缩 随机 数 产生 因子 区 间 策 


确定 当前 最 优 目标 函数 值 对 应 的 生物 体 8 pest 


采用 改进 受益 因子 策略 的 互利 共 
生 搜 索 阶段 
= 物 个 体 9;， 从 生态 3 统 中 随机 选择 生物 90,(j#i)， 两 者 进入 采用 
j ART et 生 搜 索 阶段 ， 按 式 (12)(13)(14)(19)(20) 
行 更 新 ， 产 生 新 解 9 E AAMA 9 和 8) 进行 解 的 质 
teen. 分 别 选择 较 优 解 保留 后 进入 下 一 步 。 


略 的 偏 利 共生 搜索 阶段 Y 
随机 选择 两 个 不 同 生物 9; 和 1 进入 采用 收缩 随机 数 挛 : 生 因子 区 
BERE DERE AIEE, 按 式 (21) 作 更 新 ， 产 生 新 解 Oinen 

并 与 9 进行 解 的 质量 的 比较 评价 ， 选 择 较 优 解 保留 后 进入 下 一 步 。 


寄生 阶 自 Y 
复制 生物 9,， 随 机 选取 其 让 菜 下 维度 ， 特 摘 为 搜索 人 
mus 半生 VER AIME 


S pE) Fë FE V3 ff] B 
e HENTIG 
时 es Beem 宿主 6; 


i< 


Fig.3 Flowchart of proposed imsos-based PV module models 


parameter identification 


i=i+1 


S y — 
— — ym] ea 
Æ v d8 
| 输出 光伏 组 件 模 型 的 最 优 参 数 解 向 量 2 和 对 应 的 目标 函数 值 OFewsE 
— i 
E 
Q 结束 J 
图 3 采用 所 提 ImSOS 算法 的 光伏 组 件 模型 参数 辨识 流程 


dl 各 算法 内 部 依赖 性 控制 参数 设置 
Table 1 Inner dependent parameters of various algorithms 
算法 内 部 依赖 性 控制 参数 取 值 
ImSOS 无 无 
SOS 无 无 
加 速 常数 cl 和 cz cl=c2=2 
最 小 惯性 权重 wmin 0.4 
PSO-w 最 大 惯性 权重 wmax 0.9 
最 小 速度 Vmin -10 
最 大 速度 Vmax 10 
记忆 库 取 值 概率 HMCR 0.95 
音调 微调 最 小 概率 PA Rain 0.1 
IHS 音调 微调 最 大 概率 PARmax 0.7 
微调 最 小 步 长 pwmin 1x104 
微调 最 大 步 长 pwmax 1 
对 数 螺 旋 形 状 常数 b 1 
MFO 必 敛 常数 7 [2, -1] 
随机 数 £ [r, -1] 
虫 洞 存在 最 小 概率 WEPmin 1 
MVO 虫 洞 存在 最 大 概率 WEPmax 0.2 
局 部 开发 的 准确 度 6 
SSA 系数 cl 2exp(-(411W1tmax)2) 
系数 a 2- 21t/Itwax 
对 数 螺旋 形状 常数 b 1 
WOR: 收敛 常数 7 [2,-1] 
随机 数 £ [r. -1] 
SCA 常数 a 0.8 


4.1 
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APIR 
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算 例 1 光伏 组 件 的 SDPM 模型 参数 辨识 
本 算 例 中 ,Sharp ND-R250A5 光伏 组 件 的 SDPM 模型 需 


要 辨识 五 个 未 知 参数 。 该 模型 的 五 个 未 知 参数 的 搜索 
值 分 别 为 : Tpn (A)E[0, 10], Za (AA) €[0, 1], a€[1, 2], 

E[0, 2], Ran (Q) &[0, 100]. Sz jill 
据 构成 。 根 据 实 验 工 


si 


| Ht 
Rs (Q) 
时 数据 由 36 组 电压 电流 数 
居 进 行 建 模 研 究 ， 本 文采 用 所 提 


V d 


ImSOS 算法 和 其 他 八 种 算法 来 辨识 Sharp ND-R250A5 光伏 
组 件 的 SDPM 模型 的 最 优 参 数 。 


3X2 统计 了 30 次 独立 运行 实验 中 各 辨识 


差 最 小 化 为 目标 函数 
Median, Worst 及 Std 分 别 表 示 30 次 独立 运行 实验 所 获 结果 


| 


ix. 


的 最 优 值 、 


其 中 


4 算法 以 均 方 根 误 
OFRwsE 的 结果 值 。 其 中 Best, Mean, 


平均 值 、 


中 位 值 、 最 差 值 和 标准 差 。 值 得 注意 的 
是 , 表 中 粗 体 表 示 的 数据 为 在 对 比 中 
K 2 可 以 发 现 , 对 于 Sharp ND-R250A5 光伏 组 件 的 SDPM Tz 
型 参数 辨识 结果 值 , ImSOS 算法 所 得 结 
PF 位 值 、 最 差 值 和 标准 差 均 优 于 标准 SO 算法 及 其 余 
ImSOS 4$ ik Br & 
1.13972181031247e-02, mi H. ImSOS 算法 所 得 结果 的 标准 3 


5 优 的 数据 《下 同 )。 


二 果 的 最 优 值 、 平 均值 、 
七 种 算 


x RU d UL (RON 


为 1.141686e-02， 低 于 标准 SOS 算法 及 其 他 对 比 算法 ， 这 明 


显 体现 
和 稳定 


出 所 提 ImSOS 算法 相对 于 标准 


SOS 算法 在 寻 优 精度 


t 


方面 都 有 所 提高 , 表 


明 算法 改进 具有 一 定 的 有 效 性 ， 


能 够 精确 可 靠 地 辨识 
为 了 更 直观 地 比较 ImSOS 算法 与 标准 SOS 算法 及 其 他 


七 种 人 工 智能 优化 算法 的 寻 优 


光伏 组 件 


的 SDPM 模型 参数 。 


性 能 ， 图 4 给 出 了 各 算法 得 到 


的 最 小 适应 度 值 所 对 


的 收 


取 , 各 算法 


所 提 


利 共 


部 极 值 能 力 更 强 


应 的 以 均 方 根 误差 为 目标 函数 OFRMsE 


SX iii 


线 。 从 图 中 可 以 看 出 ， 由 于 是 随机 
的 目标 初始 值 都 较 大 ， 
算法 仅 在 100 代 左 右 就 快速 收敛 到 
值 附 近 ， 并 且 在 
准 SOS 算法 及 其 余 七 种 算法 更 小 的 目标 值 。 该 结果 也 证 明了 


后 续 的 极 值 搜索 


对 初始 解 进 行 选 
随 着 迭代 次 数 增加 , ImSOS 
一 个 相对 稳定 的 目标 函数 
化 过 程 中 ， 能 搜寻 到 比 标 


ImSOS 算法 在 种 群 初始 化 过 程 中 采用 的 准 反射 学 习 机 


制 能 增加 算法 利 


ImSOS 算法 的 寻 优 能 


表 2 


FP 群 的 多 样 性 ， 提 高 了 全 局 寻 优 能 力 ， 
生 搜 索 过 程 中 采用 改进 受益 因子 策略 和 
阶段 采用 收缩 随机 数 产 生 因 子 区 间 策 略 ， 
p XR SER 
能 很 好 地 平衡 全 局 搜索 能 力 和 


算 例 


而 在 互 
在 偏 利 共生 搜索 
使 得 算法 的 跳出 局 
部 开采 能 力 。 因 此 ， 算 法 
部 开采 人 能力。 可 见 所 提 
g 较 其 他 八 种 算法 更 优 。 

寻 优 结 果 的 统计 值 


ng 


T 


zr à zz bri 


和 收敛 性 
1 中 各 算法 


Table 2 Using various algorithms 


算法 


Best 


Mean Median Worst Std 


ImSOS 1.13972181031247e-02 1.352685e-02 1.139722e-02 7.396162e-02 1.141686e-02 


SOS 1.13972181031251e-02 2.778745e-02 1.139722e-02 1.723234e-01 4.189177e-02 


PSO-w 5.229317e-02 1.180738  9.166063e-01 2.108533  Á5.454449e-0 
IHS 2.163625e-02 3.970194e-013.355801e-01 1.151087  2.669930e-0 
MFO 1.139864e-02 9.339897e-019.166063e-01 2.108533  5.780722e-0 
MVO 1.823445e-02 2.500029e-01 2.320250e-01 5.898403e-01 1.614256e-0 
SSA 1.756698e-02 2.170548e-01 1.981524e-01 4.947760e-01 1.273422e-0 
WOA 1.204897e-02 1.169997 .279503 2.108533 5.921692e-0 
SCA 7.233619e-02 3.048438e-012.266997e-01 2.108622  3.834551e-0 
表 3 给 出 了 不 同 优化 算法 进行 辨识 获得 的 最 优 参 数 如 Ion 


Isas a~ Rss Rsh 
可 以 发 现 ， 采 用 
件 的 SDPM 模型 
函数 OFrmse 值 最 小 且 
的 准确 度 , 胜 过 其 


的 值 及 相对 应 的 目标 函数 OFRMsE (Eo H 
ImSOS 算法 进行 Sharp ND-R250A5 光伏 组 
参数 辨识 时 得 到 的 最 优 参数 解 所 对 应 的 目标 


2 八 种 HHR S 


K3 


为 1.13972181031247e-02, 获得 了 最 高 
法, 按 准确 度 可 依次 排序 为 : 


ImSOS>SOS>MFO>WOA>SSA>MVO>IHS>PSO-w> SCA, 
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录用 定稿 ER od, 
这 体现 了 ImSOS 算法 增强 了 标准 SOS 算法 的 寻 优 能 力 ， 使 
其 性 能 得 到 提升 。 同 时 也 表明 采用 本 文 所 提 ImSOS 算法 可 以 
精确 地 辨识 光伏 组 件 的 SDPM 模型 参数 , 使 得 通过 所 得 最 优 


等 : 采用 改进 型 SOS 算法 的 光伏 组 件 模型 参数 辨识 


ChinaXiv 合 作 期 刊 
第 37 卷 第 4 期 


1.123868e-02， 明 显 低 于 标准 SOS 算法 及 其 他 对 比 算 法 。 特 
别 值得 注意 的 是 ，ImSOS 算法 所 得 结果 的 标准 差 为 
2.977265e-07, 远 远 低 于 标准 SOS 算法 及 其 余 七 种 对 比 算法 ， 
此 很 好 地 体现 出 所 提 ImSOS 算法 相对 于 标准 SOS. 算法 在 


参数 解 确定 的 模型 很 好 地 描述 了 实际 光伏 组 件 的 输出 特性 。 
4.5 


T T T T T T T T T 可 


目标 函数 OFrmse 


”0 200 400 600 800 1000 1200 1400 1600 1800 2000 
迭代 次 数 / 次 

图 4 算 例 1 中 各 算法 的 收敛 曲线 

Fig.4 Convergence of various algorithms for case 1 


RI 各 算法 获得 光伏 组 件 的 SDPM 模型 参数 辨识 结 


Table3 Best results of SDPM model for PV module using 
various algorithms 
算法 Is(A) Ia (uA) a R.(Q) Ra (Q) OFRMsE 
ImSOS 9.143964 1.000000 1.206865 9.854653e-03 7.759294e+01 1.13972181031247e-02 
SOS 9.143964 1.000000 1.206865 9.854653e-03 7.759290e+01  1.13972181031251e-02 


PSO-w 9.314662 3.508695e-02 1.000000 1.069221e-02 1.352126 5.229317e-02 
IHS 9.105059 6.277072e-01 1.173231 9.988120e-03 7.263959e-01 2.163625e-02 
MFO 9.142935 1.000000 1.206862 9.855259e-03 1.000000e-02 1.139864e-02 
MVO 9.233459 8.398585e-01 1.194204 9.876043e-03 3.876607 1.823445e-02 
SSA 9.166212 9.975954e-01 1.206567 9.893198e-03 9.762082e*01 1.756698e-02 
WOA 9.144157 9.237400e-01 1.201044 9.886278e-03 3.760251e*01 1.204897e-02 


SCA 9.040952 2.052154e-01 1.099011 1.007876e-02 1.972668e*01 7.233619e-02 


为 了 进一步 检验 本 文 所 提 ImSOS 算法 的 光伏 组 件 模型 
参数 辨识 质量 , 将 表 3 中 ImSOS 算法 辨识 得 到 的 最 优 参数 如 
Ip. Là. à. Rss Rs 的 值 代 入 描述 光伏 组 件 的 SDPM 模型 特 
的 方程 式 (9), 重 构 出 该 组 件 模 型 在 各 实验 测量 电压 数据 下 
对 应 的 计算 电流 数据 和 计算 功率 数据 并 绘制 出 TV、 功 率 一 
电压 (P-V) 特 性 ， 与 实验 数据 进行 比较 ， 如 图 5 所 示 。 由 图 5 
可 知 ,采用 ImSOS 算法 辨识 获得 的 最 优 参 数 重 构 出 光伏 组 件 
的 SDPM 模型 的 IV、P-V 数据 与 实验 测量 的 IV、P-V 数据 
非常 接近 ， 几 乎 难 分 彼此 , 这 表明 本 文 所 提 ImSOS 算法 获得 
的 优化 参数 具有 较 高 的 准确 性 ， 可 以 精确 有 效 地 辨识 光伏 组 
件 的 SDPM 模型 参数 。 
4.2. 算 例 2 光伏 组 件 的 DDPM 模型 参数 辨识 

为 进一步 验证 本 文 所 提 ImSOS 算法 辨识 光伏 组 件 模 型 
参数 的 有 效 性 和 优越 性 ， 本 算 例 中 ， 模 型 需 辨 识 的 未 知 参 数 
从 算 例 1 的 五 个 增加 至 七 个 , 其 辨识 求解 难度 更 大 , B Sharp 
ND-R250A5 光伏 组 件 的 DDPM 模型 需要 辨识 七 个 未 知 参 数 。 
该 模型 的 七 个 未 知 参 数 的 搜索 范围 取 值 分 别 为 : Iph (A) € [0， 
10], Ka (4A) €[0, 1], Zo (4A) El[O, 1], aie [1, 2], a2 e [1, 2], 
Rs (Q) € [0, 2], Rsa (Q) € [0, 100]. 实验 测量 I-V 数据 同 算 例 1。 
同样 采用 本 文 所 提 ImSOS 算法 和 上 述 8 种 算法 来 进行 Sharp 
ND-R250A5 光伏 组 件 的 DDPM 模型 最 优 参数 辨识 研究 。 
在 对 各 辨识 算法 以 均 方 根 误差 最 小 化 为 目标 函数 
OFrMsE 独立 运行 30 次 后 ,并 对 其 结果 值 分 别 进行 统计 分 析 ， 
如 表 4 所 示 。 
表 4 可 以 清楚 地 看 出 , 对 于 Sharp ND-R250A5 光伏 组 
件 的 DDPM 模型 参数 辨识 结果 值 ，ImSOS 算法 获得 的 结果 
的 最 优 值 、 平 均值 、 中 位 值 、 最 差 值 和 标准 差 均 优 于 标准 SOS 
算法 及 其 余 七 种 算法 。 其 中 Imsos 算法 所 得 结果 的 最 优 值 为 


HE 


稳定 寻 优 能 力 和 准确 性 方面 都 有 所 提升 ， 表 明 算 法 改进 具有 
一 定 的 有 效 性 , 在 辨识 光伏 组 件 的 DDPM 模型 最 优 参数 时 具 
有 更 为 可 靠 的 准确 性 和 稳定 性 。 
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(b) P-V 特性 

图 5 光伏 组 件 的 实验 测量 数据 与 ImSOS 算法 辨识 其 SDPM 模 

型 的 计算 数据 对 比 


Fig. 5 Comparisons between measured data and calculated data 


by proposed imsos for SDPM model 

图 6 直观 地 给 出 了 各 算法 寻 优 得 到 的 最 小 适应 度 值 所 对 
应 的 以 均 方 根 误差 为 目标 函数 OFrmse 的 收敛 曲线 。 从 图 6 
中 可 以 发 现 ， 各 算法 的 目标 初始 值 较 大 ， 随 着 迭代 次 数 的 增 
加 ，ImSOS 算法 仅 在 200 代 左 右 就 快速 收敛 到 一 个 相对 稳定 
的 目标 函数 值 附 近 ， 并 且 在 后 续 的 极 值 搜索 优化 过 程 中 ， 能 
搜寻 到 比 标准 SOS 算法 及 其 他 算法 更 小 的 目标 值 。 该 结果 
证 明了 所 提 ImSOS 算法 改善 了 标准 SOS 算法 的 寻 优 性 能 ， 
能 很 好 地 平衡 全 局 探索 能 力 和 局 部 开采 能 力 。 可 见 所 提 
ImSOS 算法 的 寻 优 能 力 和 收敛 性 能 较 其 他 八 种 算法 更 优 。 

表 5 清楚 地 给 出 了 各 优化 算法 进行 辨识 获得 的 最 优 参 数 
如 fon. fus dem. di. a2. Rs. Ron 的 值 及 相对 应 的 目标 函数 
OFRwsE 值 。 由 表 5 可 知 ， 采 用 ImSOS 算法 进行 Sharp 
ND-R250A5 光伏 组 件 的 DDPM 模型 参数 辨识 时 得 到 的 最 优 
参数 解 所 对 应 的 目标 函数 OF nus 值 最 小 且 为 1.123868e-02， 
获得 了 最 高 的 准确 度 ， 明 显 优 于 其 余 八 种 辨识 算法 ， 按 准确 
度 可 依次 排序 为 : ImSOS>SOS>MFO>WOA>SSA>MVO> 
IHS>PSO-w>SCA， 这 充分 体现 了 ImSOS 算法 改善 了 标准 
SOS 算法 的 寻 优 能 力 ,特别 在 本 算 例 中 求解 问题 维 数 增加 时 ， 
其 性 能 的 提升 体现 得 更 加 明显 。 同 时 也 表明 采用 本 文 所 提 
ImSOS 算法 在 精确 地 辨识 

光伏 组 件 的 DDPM 模型 参数 时 具有 一 定 的 优越 性 , 使 得 
通过 所 得 最 优 参数 解 重 构 获 得 的 模型 与 实际 光伏 组 件 输出 特 
性 很 好 地 吻合 。 
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录用 定稿 Ax, 


dA 算 例 2 中 各 算法 寻 优 结果 的 统计 值 


Table 4 Using various algorithms 


算法 Best 
ImSOS 1.123868e-02 


Mean Median Worst Std 


1.123881e-02 1.123868e-02 1.123986e-02 2.97726Se-07 


等 : 采用 改进 型 SOS 算法 的 光伏 组 件 模 型 参数 辨识 


till, 


利 共生 搜索 阶段 采用 了 改进 受益 因子 策略 以 使 每 个 生物 在 它 
门 相互 作用 中 获得 相等 的 部 分 益处 ， 提 高 了 算法 的 稳定 性 ; 
在 偏 利 共 生 搜 索 阶 段 采 用 了 收缩 随机 数 产 生 因 子 区 间 策 略 ， 
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提高 了 算法 的 求解 精度 和 加 快 了 算法 的 收敛 速度 ; EH. 


SOS 1.123895e-02 4.435637e-02 1.132074e-02 8.808801e-01 1.592047e-0 把 慢 索 范围 控制 在 合理 区 间 ， 加 快 了 算法 的 收敛 速度 的 同时 
PSO-w 2.314774e-02 1.010029 9.157699e-01 2.108533 4.090412e-0 岂 保 证 了 其 求解 精度 。 其 次 ， 给 出 TH} 所 提 ImSOS 算法 求 
aro ier cima rca ro mca PATTERN LV BARMERA SAONA 
MVO 1.235726e-02 1.115238e-01 7.324502e-02 3.716043e-01 1.070155e-0 本 步 又 及 实现 流程 ,最 后 ,利用 所 是 ImSOS 算法 对 实际 Sharp 
SSA  1.207423e-02 5.647215e-02 4.726401e-02 1.915108e-01 3.946225e-02 ND-R250A5 光伏 组 件 进行 J 实验 验证 。 通过 与 标准 SOS, 
WOA 1.168769e-02 1.453825 1.277271 6.671749 1.888110 PSO-w, IHS, MFO, MVO, SSA, WOA 和 SCA 8 FIUBUSUS 
SCA  4.543367e-02 1.719073e-01 1.509657e-01 9.321523e-01 1.514787e-0 能 优化 算法 进行 对 比分 析 , 结果 证 明了 所 提 ImSOS 算法 增强 
se 了 标准 SOS 算法 的 寻 优 能 力 , 使 其 性 能 得 到 了 提升 。 同 时 也 
: — sos] | 表明 了 所 提 ImSOS 算法 在 光伏 组 件 的 SDPM 模型 和 DDPM 
" Z Pow | | 模型 参数 辨识 的 有 效 性 和 优越 性 ， 为 准确 可 靠 地 辨识 光伏 组 
d. — MP | 件 模型 参数 提供 种 新 的 有 效 方法 ， 以 期 为 实际 中 光伏 组 T 
E. —s || 2 UU 性 能 评估 以 及 故障 检测 等 方面 研究 提供 有 力 的 
性 — — SCA 村 。 
=>? s | 然而 本 文 研 究 并 未 考虑 不 同 遮 戎 等 外 部 因素 情况 下 光伏 
"i ph ] 组 件 遮挡 模型 的 建立 以 及 模型 参数 的 精确 辨识 ， 后 期 研究 将 
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和 迭代 次 数 /次 

图 6 算 例 2 中 各 算法 的 收敛 曲线 

Fig.6 Convergence of various algorithms for case 2 

为 了 更 形象 地 说 明 本 文 所 提 ImSOS 算法 的 光伏 组 件 模 

型 参数 辨识 的 准确 性 , 将 表 5 中 ImSOS 算法 辨识 得 到 的 最 优 

参数 如 Dos. Fai. Ja. ai. azs Rs. Ran 的 值 代入 描述 光伏 组 

件 的 DDPM 模型 特性 的 方程 式 (10)， 重 构 出 该 组 件 模 型 在 各 
实验 测量 电压 数据 下 对 应 的 计算 电流 数据 和 计算 功率 数据 

绘制 出 I-V、P-V 特性 ， 与 实验 数据 进行 比较 ， 如 图 7 所 示 。 

由 图 7 可 以 直观 地 看 出 ,采用 ImSOS 算法 辨识 的 最 优 参数 习 

构 出 光伏 组 件 的 DDPM 模型 的 ITV、P-V 数据 与 实验 测量 的 

I-V, P-V 数据 同样 非常 接近 ， 有 具有 很 好 的 吻合 度 ， 这 表明 本 

文 所 提 ImSOS 算法 获得 的 优化 参数 具有 较 高 的 精度 , 尽管 该 

算 例 的 组 件 模 型 中 未 知 参 数 增多 了 ， 辨 识 难度 增 大 了 ， 同 样 
可 以 准确 有 效 地 辨识 出 光伏 组 件 的 DDPM 模型 最 优 参数 。 

表 5 各 算法 获得 光伏 组 件 的 DDPM 模型 参数 辨识 结果 
Table 5 Best results of DDPM model for PV module using various 
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Tp 


algorithms 


算法 Im(A) [a (pA) (pA) ai a R: (9) Ra (Q) OFRwsE 


会 重点 考虑 这 方面 的 问题 。 此 外 ， 下 一 步 也 可 继续 探索 将 本 


文 所 提 ImSOS 算法 


E 广 到 其 他 实际 工程 优化 问题 求解 的 可 


ImSOS 9.145875 9.212180e-02 1.000000 1.213483 1.213484 9.824284e-03 1.000000e+02 1.123868e-02 


SOS 9.145883 7.874803e-01 3.053844e-01 1.211853 1.218075 9.824266e-03 9.999999e+01 1.123895e-02 


PSO-w 9.209539 1.000000 — 2.978648e-01 2.000000 1.123047 1.018894e-02 3.388706 . 2.314774e-02 


IHS 9.164499 7.092563e-01 1.000000 1.191551 1.384045 9.855943e-03 9.476771e*01 1.619519e-02 


MFO 9.144779 1.000000 1.000000 — 1.832347 1.207152 9.844784e-03 1.000000e-02 1.133178e-02 
MVO 9.150307 9.354702e-01 2.792301e-04 1.201868 1.999662 9.888661e-03 5.756872e*01 1.235726e-02 
SSA 9.144893 8.465336e-01 3.182345e-01 1.195376 1.526311 9.883436e-03 2.883042e+01 1.207423e-02 
WOA 9.139977 9.667885e-01 9.925868e-01 1.204515 1.911408 9.848554e-03 9.924384e*01 1.168769e-02 


SCA 9.093958 1.795627e-24 6.331248e-01 1.360216 1.173024 9.802448e-03 5.319606e*01 4.543367e-02 


5 ”结束 语 


本 文 提 出 了 一 种 采用 改进 型 SOS 算法 的 光伏 组 件 模 型 
参数 辨识 方法 。 首 先 ， 为 提高 标准 SOS 算法 的 寻 优 能 力 ， 提 
出 了 新 的 改进 型 共生 生物 搜索 算法 , 记 做 ImSOS 算法 。 该 算 
法 采用 了 三 个 主要 策略 来 改善 标准 SOS 算法 的 优化 性 能 , 很 
好 地 平衡 了 算法 的 全 局 探索 能 力 与 局 部 开发 能 力 。 该 算法 在 
标准 SOS 算法 的 生物 种 群 初始 化 阶段 采用 了 准 反射 学 习 机 
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Fig.7 Comparisons between measured data and calculated data by 


proposed imsos for DDPM model 
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